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Abstract. In this paper we will discuss analytically the perturbations created from a slowly 
rolling subdominant spectator field which decays much before the end of inflation. The 
quantum fluctuations of such a spectator field can seed perturbations on very large scales 
and explain the temperature anisotropy in the cosmic microwave background radiation with 
moderate non-Gaussianity, provided the relevant modes leave the Hubble patch while the 
spectator is slowly rolling. Furthermore, the perturbations are purely adiabatic since the 
inflaton decay dominates and creates all the Standard Model degrees of freedom. We will 
provide two examples for the spectator field potential, one with a step function profile, 
and the other with an inflection point. In both the cases we will compute higher order 
curvature perturbations, i.e. local bispectrum and trispectrum, which can be constrained by 
the forthcoming Planck data. 
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1 Introduction 

Inflation is the most well known paradigm for the early universe which has confronted the 
observations very well [1]. It is a dynamical mechanism which stretches the quantum per- 
turbations on very large scales, while making the universe big and diluting all matter [2-4]. 
For a review on inflation models, see [5]. During inflation there could be more than one field 
participating in the dynamics - they may drive inflation collectively as in the case of assisted 
inflation [6], or there could be one or many subdominant fields who do not directly influence 
the overall dynamics. Such a field might even decay early on while inflation is still going on, 
and we may regard it as a spectator field. In principle there could be many spectator fields . 

If the potential for a spectator field is sufficiently flat before it decays, it can slow roll, 
and during this phase it can accumulate Hubble induced quantum fluctuations. These fluctu- 
ations can be imprinted in the temperature anisotropy of the cosmic microwave background 
(CMB) radiation. In the simplest scenario we may assume that the inflaton fluctuations are 
subdominant, and a single spectator field is present in the universe besides the inflaton which 
drives inflation. 

The role of a subdominant field during inflation has been discussed before in the context 
of a curvaton scenario, where a light field seeds the fluctuations and then decays after inflation 
comes to an end [10-13]. Very recently we have discussed a novel possibility of creating 

1 The spectator mechanism is different from the multi-field inflation scenarios discussed in [7-9]. 
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perturbations from a spectator field which decays while inflation is still going on [14]. There 
are some crucial differences between these two scenarios which we highlight below. 

Since inflation dilutes all matter, it is then necessary to create all the Standard Model 
degrees of freedom after the end of inflation, for a review on reheating, see [15]. In the 
context of a curvaton scenario one possibility is that the curvaton produces all the relevant 
matter and dominates the energy density of the universe while decaying [16-18]. This is 
however difficult in most realistic scenarios, where we also require the curvaton to decay 
solely into the Standard Model degrees of freedom. Another possibility is that the curvaton 
subdominates the universe, so the relic isocurvature perturbations is weaker than the 5% 
observational bound [1]. Otherwise the decay products of the inflaton's and the curvaton's 
have to thermalize, but this is a daunting task. Usually implicit assumptions are being made 
that the decay products would thermalize, but there are few setups which can explicitly 
realize this in reality [19]. A subdominant curvaton also leads to a significant enhancement 
in non-Gaussianity [12, 13]. 

On the other hand if a field decays much before the end of inflation, it will never influence 
the thermal history of the universe and will become a spectator field. Its decay products will 
be redshifted away during inflation and the inflaton will be solely responsible for creating all 
the matter in the universe. Meanwhile such a spectator field could still be responsible for 
seeding the CMB anisotropy as discussed first in [14], provided the relevant scales for the 
CMB leave the Hubble patch before the spectator field decays. In this respect there will be 
only adiabatic perturbations and no isocurvature perturbations in the universe. 

In order for this scenario to work, obviously it is important that the spectator field's 
potential is very flat during inflation before it decays. There are no dearth of such fields, 
their origin could come from anywhere, i.e. from the visible sector (beyond the Standard 
Model sector whose fields carry the Standard Model charges), or from a hidden sector (fields 
which do not carry the Standard Model charges). Such a spectator does not even have to 
couple to the Standard Model degrees of freedom. All the onus will be now on the inflaton's 
coupling to the Standard Model degrees of freedom for creating the right thermal history of 
the universe. In this case inflation could be solely driven within the visible sector [20, 21], or 
it may arise from any hidden sector with a specific coupling to the matter fields [22-24]. 

To be specific, we will derive general formulae for a generic potential of the inflaton and 
the spectator field, with the help of 5N formalism developed in [25, 26], for a pedagogical 
discussion, see [27]. Then we will consider an example where the spectator field has a sharp 
edge in the potential, where the slow roll conditions are naturally violated, and the final 
expressions are greatly simplified. In the second example we will consider a smooth transition 
from the slow roll phase to the oscillatory/decay phase of the spectator field. This can happen 
in the case where the spectator potential has an inflection point. The inflection point is where 
the second derivative of the potential vanishes but not the first or third derivatives. This 
gives a gentle slope for the spectator field to roll down the potential before the oscillation or 
the decay begins. 

Before the spectator field decays, there exists an entropy perturbations due to the fluc- 
tuations in the spectator and in the inflaton field, although the inflaton's perturbations are 
subdominant. Both the examples provide mild or large local bispectrum. In the first exam- 
ple, the spectator has a flat and smooth potential, so the mild bispectrum is mostly sourced 
by the conversion from the short lived entropy perturbations to the curvature perturbations, 
which becomes non-Gaussian after the spectator decays. The second example has a curved 
potential for the spectator field, due to the inflection point, and therefore can generate signif- 
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icant local bispectrum from the potential curvature and from the conversion of the entropy 
perturbations into the curvature perturbations. In both cases, the non-Gaussianity is sourced 
by the entropy perturbations as it was highlighted in the context of non-perturbative decay 
of the inflaton [28, 29]. Eventually, the spectator field or its decay products are redshifted 
away during inflation, leaving no trace of isocurvature perturbations. In both the examples, 
the cosmological observables are predicted within our current observational limits. Note that 
this scenario generates non-Gaussianity in a different way than the multi-field inflation mod- 
els, where large non-Gaussianities are expected due to non-trivial end-of-inflation boundary 
conditions, for example, see [30]. 

We also have to keep in mind that at the Hubble exit of the CMB's relevant scales, the 
spectator must have its classical slow roll dominating over its quantum fluctuations. There- 
fore the potential for the spectator field cannot be extremely flat or extremely subdominant. 
Otherwise this would give rise to a stochastic phase for the spectator field [31], which we 
will not explore here. This condition further puts a model independent lower bound on the 
energy density ratio between the spectator and the innaton when the spectator decays. This 
lower bound will automatically prevent the generation of a very large local bispectrum or 
trispectrum. 

In Section 2, we will discuss the general setup. In Section 3, we will calculate the 
curvature perturbations at the linear and the higher orders. In Section 4 and Section 5, we 
will provide examples for the spectator field with a step potential and an inflection point 
potential, respectively. We will conclude this paper in Section 6. A list of variables are 
summarized in Section A. 

2 General setup 

Let us consider a generic potential with the inflaton, <j), and the spectator, a. They don't have 
any interactions except their minimal couplings to gravity, so their potential has a general 
form: 



Here we assume both the fields are canonical scalar fields, which always satisfy V(<p) 3> U (a). 
We also assume the spectator field a ends slow roll well before the end of inflation but after 
the Hubble exit of the pivot scale. This will give rise to two subsequent phases for a during 
inflation as shown in Fig 1: 

1. Phase I: The inflaton <fi leads inflation. Both eft and a are slowly rolling. This phase 
ends as the second order slow roll condition for a breaks down. We assume that the 
relevant perturbations for the CMB are leaving the Hubble patch in this phase 2 . 

2. Phase II: When a ends slow roll (depicted here by point "c"), the inflaton <p still 
dominates inflation under slow roll. Then a either oscillates around the minimum of 
it's potential, or decays instantly once its slow roll has been terminated. In either case, 
a or its decay products are being redshifted away swiftly during this phase, and can be 
regarded as a perfect fluid with a constant equation of state w. For this reason, several 
e- folds after the beginning of this phase, the dynamics reduces to that of a single field 
inflation. 

2 The observed pivot scale actually has a width of several e-folds. Here we consider every mode separately 
so the window is not shown in Fig 1. 




(2.1) 
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* Phase I c Phase II e 

Hubble exit cr ends slow roll End of inflation 

Figure 1. A schematic timeline of the universe's evolution is shown above. The solid lines are the 
background evolutions of the energy densities of the inflaton, the spectator and the total contributions. 
The green dashed lines are the total energy densities of the universes with perturbed a . The two phases 
of a evolution are separated by the phase boundary when a breaks the slow roll condition. The "*" 
denotes the epoch when a fluctuations leave the Hubble patch, "c" denotes the boundary between 
the two phases, and "e" denotes the end of inflation. 



We assume that the primary curvature perturbations are generated from the a field and 
the inflaton's perturbations are negligible. During phase I, the energy density of a contributes 
to the Hubble rate, albeit small. But in phase II, its energy density or its decay product's 
contribution is quickly redshifted away by the ongoing inflationary expansion. Perturbations 
of the a field which are depicted in Fig 1 by the dashed lines induce fluctuations of the phase 
boundary itself, which converts the a perturbations into curvature perturbations. 

The cosmological observables can then be calculated using the 8N formalism [25, 26]. 
According to 5N formalism we can regard separate universe patches as being perturbed as a 
whole by the perturbation modes at the Hubble exit. The different initial field perturbations 
(5(7* at different universe patches then lead to different 5N, the perturbations in the number 
of remaining e-folds of inflation, which directly link to the amount of curvature perturbations 
as defined by £ [32] , for a review see [33] . 

For modes leaving the Hubble patch during phase I, the fluctuations, 5a*, of a flat 
enough spectator field are Gaussian random perturbations with the power spectrum 

= Hi (2-2) 

where H is the Hubble expansion rate of inflation, and "*" indicates the Hubble exit corre- 
sponding to the pivot scale. Imagining N as a function of cr*, we can write its perturbation 
as 

SN = N a 5o* + \N aa (5a*) 2 + \N a(J(7 {5a*f + 0{{5a*f), (2.3) 
2 o 
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where N a = dN/da*, N acr = d 2 N/da* and N a(Ta = d 3 N/da*. Therefore the power spectrum 
for the curvature perturbations can be expressed at the first order as 

P ( = P 5N = N 2 P 5a ,. (2.4) 

The deviation from a Gaussian distribution of curvature perturbations is expressed 
order by order. The two lowest order parameters, namely the local bi-spectrum (/nl) and 
the tri-spectrum (<7nl)> have their relative local strengths characterized by 3 

, 5AU 25 N aaa 

= 6 ~N%~ ' 5NL = 54lvr- (2 - 5) 

Therefore in order to obtain the observational predictions, we need to work out the derivatives 
of N as in Eq. (2.3). In accordance with the 5N formalism, we will use dN = —Hdt. 

3 Curvature perturbations 

In order to understand the effects of the perturbation 5a*, we will consider the two phases 
separately and study how their numbers of e-folds are perturbed. After finishing the calcu- 
lation at first order, we will use the results to obtain answers for the higher orders. Note 
that in our perturbative calculation the only initial perturbation is 5a*, while all the other 
5's are the induced perturbations by the initial perturbation 5a*. 

3.1 Perturbations from phase I 

We start from the background equations of motion for both the fields with the total potential 
Eq. (2.1). Under the slow roll conditions, we use N as the remaining number of e-folds of 
inflation as a proper time, and rewrite their equations of motion as 4 

dl MlV da_ Mp* 

dN 8vr(f/ + y)' dN 8tt(U + V)' [ ' ' 

where primes on the potentials denote derivatives w.r.t the corresponding fields, i.e. V = 
dV/d(j) and U' = dU/da. These indicate a simple relation 

# _ da 

V ~ XT' 1 j 

whose integrated form is 

%-£% 

In the presence of an initial field perturbation 5a* , Eq. (3.3) will also be perturbed, 
which leads to the perturbations of 5(f) and 5a at any time, with the relation 

_ _ 5a*^ _ 5a_ 
V U*' u r { ' ' 

3 The other parameter for the strength of trispectrum, tnl, is reduced to be proportional to f^ L , when 
the quantum fluctuations of one field dominate over the rest. 

4 Our formalism for phase I is similar to Ref [34]. However this formalism cannot be applied to phase II 
directly. 
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The phase boundary between phase I and phase II is the breakdown of the second order 
slow roll condition for a 

f = M^U" + 8Tr(U + V), (3.5) 

so for the case we consider, at the phase boundary "c" we will have f c = 0. Since the phase 
boundary relation remains the same regardless of how the separate universe is perturbed, in 
the perturbed universe we would have the field values at the phase boundary also perturbed 
by the amount 5<p c and 5a c , so it remains on the phase boundary 

s f\ c = UcHc + facbVc = 0, (3.6) 

where 

M*U C '" + torU e '. (3.7) 



^V c , f ac = — 



Taking Eq. (3.4) at the phase boundary "c", together with Eq. (3.6), we are able to derive 
the induced field perturbations on the phase boundary by the initial perturbation 5cr* 

6(t> c = -^(1-8)5^, Sa c =^95a„ (3.8) 

where 
Here 

ia = 64^{U + W (3 ' 10) 
is the third order slow roll parameter for a and is typically of order 1 or bigger for the case 
we consider. It can even be much greater than unity near the sharp edges of the spectator 
potential. On the other hand, the first order slow roll parameters for both the fields are 

M 2y>2 M 2 V ,2 

^ = iWTVF « ^ s itov + vy < L (3 - n) 

Starting from "*", the Hubble exit of the mode we are interested in, the number of 
e-folds in phase I can be written as an integrated form on the uniform <fi slicing 

The perturbation in iVj then becomes 

On uniform </> slicing we always have 5(f) = 0, so Eq. (3.4) becomes 5a = U'5a*/Uj . This 
simplifies the integral in Eq. (3.13) to 

Consequently, we get the final expression for the perturbed number of e-folds in phase I as 

r at 8n(U.-U c + (l-e)(U c + V c )) x 

5N, = j^-, 6*.. (3.15) 
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3.2 Perturbations from phase II 

As explained in Section 2, during phase II we have a single field inflation with an additional 
perfect fluid component whose energy density came from either the decay products or the 
oscillations of a. The Hubble expansion rate is given by 



3Mf 



(U ce S^)(N- Nc) + = ^M^ (3 . 16) 



where the second half comes from the slow roll equation of motion for 0. Here U c means the 
value of U (a) at the phase boundary c, and w is the equation of state for the perfect fluid. 

It is obvious that Eq. (3.16) is actually a first order differential equation between N and 
<j). Its exact solution is given by 

N2 = N c — N e = n((f> c , 4> e ) + j r ln^^, (3.17) 

3(1 + w) 1 — r 

where 

Z^ 1 8irV 

n(<^ 2 )=/ — (3.18) 

is the number of e- folds of inflation when eft serves as the only component of the universe, 
while the second term is the contribution from the perfect fluid. In Eq. (3.17) we have defined 

as the energy density ratio of the perfect fluid at the phase boundary, and 

24(1 + w)irV c f*° e -3(l+™M<Pc,4>) 
■ = 1 - -^-J^- 1 m # « 1 (3.20) 

will be calculated in Section B. Note that a is order of the slow roll parameter e^ c <C 1, and 
hence a <C 1. 

In a perturbed universe with Jcr*, field perturbations 8<f) c and 5o~ c will be generated on 
the phase boundary. They will give rise to the perturbations of the above parameters 

HtcA) = ^fjSfa, (3.21) 
5r = 1 ((1 - r)U c '5a c - rV c '5<j) c ) , (3.22) 

U c "T ' c 

and 

5a = -\Jl-a)— + a J^fiyl J 6 ^ ( 3 - 23 ) 

Feeding them back into Eq. (3.17) gives the perturbation in the number of e- folds for phase 
" A at 1 ( 8W C (l-a)U c ' , \ r „ 9A , 

^ 2 = T^Tr {mIV7 5(Pc + m + w )(u c + v c ) 6(Tc ) ■ (3 ' 24) 
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3.3 Total perturbations 

Adding up Eq. (3.15) and Eq. (3.24) gives the total perturbation in the number of e-folds 
generated by the initial perturbation Sa* 

5N = 5Ni + 5N 2 = N a 5a*, (3.25) 

where 

l-a f 8nU, / U c ' 2 8ttU c \ \ q(l - r) 8tt(U, - U c ) 

a l-ar\M 2 UJ \3(l + w)(U c + V c )U*> M*Uj)) 1 - or M 2 UJ ' 

(3.26) 

When the potential U(a) is sharp enough at the phase boundary, the conditions S> e^ c 
and r^o-c 3> t^, c t ac are satisfied. This is the case we will be analyzing in this paper. Under 
these conditions, Eq. (3.26) is dominated by the very first term. It will simplify to 

We see that Eq. (3.27) comes in agreement with the results obtained in Ref. [14]. When the 
spectator dominates the curvature perturbations, we need the spectator's potential should 
be relatively flatter than that of the inflaton's 5 , i.e. [/*/{/*' S> V*/VV . 

The above calculations would then give the power spectrum of the curvature perturba- 
tions 

WH 2 U 2 H 2 s 2 r 2 

where 

•=£±£<1 (3-29) 

is the ratio between the total energy densities at the phase boundary and at the Hubble exit. 
The spectral tilt is given by 6 

<91nP c ±s a £ a * 

n s - 1 = <vat" = ~ 2e 0* + 2 Vv* , (3.30) 

oJS sr 

where 

Sa = jy- < 1 (3.31) 

is the energy density ratio of a between the phase boundary and the Hubble exit. 

Note from Eq. (3.28), we can obtain a constraint on e a * from the observed P^, leading 

to 

_ H 2 s 2 r 2 

6(T * ~ 7TS 2 a M 2 P ( - (3 ' 32) 

By plugging Eq. (3.32) back into Eq. (3.30), we see that for the inflation energy scale < 
10 _5 Mp, the last term in Eq. (3.30) is always negligible compared to the observed spectral 
index, n s ~ 0.96 [1]. So we are left with 

n s - 1 = -2e0* + 27? CT *. (3.33) 



This is because when the inflaton totally dominates the curvature perturbations, we have a similar 
expression 

6 The — 2eo- „ term is neglected because it is much smaller than the last term in Eq. (3.30). 
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With the help of Eq. (3.30), the running of spectral tilt can be shown as 
dn s dn s 1 



dink dN 2 

where 



(n s - l) 2 - 2^3e 2 - 2(e^*^* + e a *rj a *) - t? 2 ^ + (3.34) 

e =|^ = 6, + e<7 (3.35) 

is the total slow roll parameter for inflation, and e* ~ e^* can be taken in Eq. (3.34) for 
< KT 5 M p . 

3.4 Higher order perturbations 

Higher order curvature perturbations, i.e. the local bispectrum, /nLj and the trispectrum, 
require the expression of a, which is derived in Section B. 

The strength of the local bispectrum /nl can be derived from taking the derivative 
<9/<9(7* on Eq. (3.26). The leading terms are 

" 6 N 2 ~ 3r V 3(1 + w) & \ UJ U 
+ 2Os 2 r e CTC £0 C / _ 4 + 3w _ X^A + 5sle a * 5s a r] a 



9(1 + w)r 2 ^ ac \ £ ac i 2 ac J 3s 2 r 2 6sr 
+ higher order terms, (3.36) 

where the third and fourth order slow roll parameters £ CT and X a are defined as 8 

M*U'U"' _ M^U' 2 U"" 

^ = (8ir(U + V)) 2 ' K = JMU + VW' (3 ' 37) 

The third order derivative N acra = d 2 N a /da 2 can be calculated in the same way and 
also #nl- According to Eq. (2.5), we obtain the leading order trispectrum of curvature 
perturbations as 

(3.38) 



where 



A = €<t>c + ^ 2r^ c H 3e0 C 2770c + 



S / V S J V C<TC 6, 

2 / 2 3(l+™ + A CTC ) 3(4 + 3™)A <JC -x <TC , 3A 2 



2 



><TC 



+ 24 c (6--- - ■~-"^ + v — + (3.39) 

SO"C Sere Sere 



and 



M*U' 3 U""' 

*> s (RcTTVF (3 ' 40) 



7 We take series expansion for the parameters e^* -C 1, e^c <C 1, <C 1, |f?<£ c | <C 1, e a * <C 1, fac <S 1, 
r?tT*| <C 1, r <C 1 and by writing down only the leading order contributions for the expressions of /nl and 
(/nl • The slow roll parameter does not have to be much smaller than 1, but in most cases it has the order 
of e| c . Therefore here we also take it as a small quantity. 

8 Note that the third and higher orders slow roll parameters of a (e.g. £ CT and A CT ) typically vary from O(l) 
to (positive or negative) infinity depending on the actual potential (7(a). Therefore we do not take them as 
infinitesimals and do not take series expansion on them in Eq. (3.36) as well as in Eq. (3.38). 
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is the fifth order slow roll parameter for a. 

Here we pause and discuss the validity of the perturbation theory itself. In all the 
above calculations, we have assumed that the classical slow roll dominates over the quantum 
fluctuations of a for the relevant scales. Here we focus on comparing the strengths of quantum 
fluctuations and classical slow roll, ensuring the validity of the above calculations. 

A dominant classical slow roll requires the field displacement of a to be larger than the 
quantum fluctuations in one Hubble time, i.e. H*/2ir < da/H*dt, or 

P fo .<(^) 2 . (3.41, 

Multiplying both sides by N%, we convert the l.h.s into the power spectrum of the curvature 
perturbations P^. A simple exercise then yields a model independent lower bound on r 

r>—JP~r. (3.42) 

8 

Since from observations P^ ~ 2.5 x 10~ 9 , we obtain a lower bound r > 5 x 10~ 4 from a simple 
estimation by using s a = U c /U* 1 and s = (U c + V C )/(U* + V*) ~ 0.1. 

From Eq. (3.36) and Eq. (3.38), we see that /nl and #nl depends very much on the 
magnitude of r. Therefore a lower bound on r typically becomes an upper bound for /nl 
and <?nl- in many cases this constraint automatically prevents the possibility of a very large 
non-Gaussianity. 



4 A step potential for the spectator field 
4.1 General predictions 

Since the spectator field decays during inflation and all the matter are created by the inflaton, 
we illustrate a very simple potential for U(a), which is given by 

U{a) = Ste P (a-a Ha) = { U ^ ^ > ^ (4.1) 

where we assume u(o~) is flat and smooth enough to accommodate slow roll for a > ctq. The 
step function is a limiting case of the hyperbolic tangent function 

1 + tanh kx ... 
Step(x) = lim . (4.2) 

k— >+oc 2 

We also assume that in the early universe a would be initially displaced at a > do, 
slowly rolling down the potential u{a). The field a ends slow roll at a c ~ do, because of 
the violation of the second order slow roll condition. At this point its energy density is 
transferred to the ideal fluid with a constant equation of state w. On the other hand when 
a ends slow roll, it still remains on the u(o~) potential which ensures tanh k(a c — do) ~ 1, i.e. 
k{a c — o"o) 3> 1. Therefore the derivatives of the tanh function at the point a c has a leading 
contribution 

dntanh d k ^- ao) « (- 2 y^k n e- 2k ^ o), for n = 1, 2, . . . under k(a - o ) » 1. (4.3) 
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With the phase boundary condition / = from Eq. (3.5), we find the deviation of a c 
from o"o as 

1 k 2 M 2 U c 

ac - ao = Tk ln Mu7W7Y (44) 

The higher order slow roll parameters then become 

i ac = 2kM p <^, (4.5) 

V 7T 

\ ac = -2k 2 M 2 ^ , (4.6) 

^ c = 2^M p 3 (v) f - (47) 

For a step function potential whose edge is infinitely sharp (i.e. k —> oo), and the 
potential u(a) is very flat and smooth (i.e. e a * *C 1, r]a* "C 1, and s CT » 1), the local 
bispectrum Eq. (3.36) is dominated by the first term, where we also have £ ac —> oo. The 
bispectrum and the trispectrum are therefore given by 

10e^ c (2e^ c - r^ c ) 

JNL ~ TTfr- s , (4.8) 

9 1+ io)r 



25gg* _50e^c^c 
54s 2 r 2 81(1 + iu)r 2 



9NL ~ "77^ + owi i 1^2 " ( 4 - 9 ) 



From Eq. (4.8) we find an agreement with [14]. 
4.2 Double quadratic potentials 

As a naive example, we consider both fields to have quadratic potentials. The potential U(a) 
also has a step function for a sharp transition during inflation. In this case, the potentials 
are written as 9 

V(</>) = m 2 .^ 2 , U[a) = m 2 a a 2 Step(cr - a ). (4.10) 

Inflation is driven by V(4>). In the beginning when the relevant perturbations leave the 
Hubble patch, we have cr* > ctq, so it stays on the plateau and rolls down slowly. When the 
spectator field reaches a c = ao, the sudden change in the potential would break the second 
order slow roll condition for a, ends its slow roll, and let it decay instantly to radiation as 
we assume. The radiation is quickly diluted away by inflation. 

For this model, we have a total of 4 free parameters — m^, m a , ao, and N c - the number 
of e-folds from the phase boundary "c" to the end of inflation. The overall energy scale only 
affects the power spectrum of curvature perturbation by 

,2^2 



16(2iV* + l)m 2 y Q 



c 3M 4 v ' 



Therefore we can fix to the observed value and hence reduce the number of parameters to 
3. We want a to dominate the curvature perturbations, which requires ao to be large enough 
((To 3> (/)*)■ However as long as this condition is satisfied, the value of ao hardly changes 
the cosmological predictions of the model. We are left with only two free parameters which 



9 Because a rolls very slowly during phase I, we typically expect <r» — ao <C ao- This means in phase I the 
effective potential for a is actually U(a) — m^ao + 2m^a (a — ao). 
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Log 10 r 

Figure 2. The local bispectrum strength /nl is shown for a spectator held with a step potential. 
The yellow shaded region is excluded in the parameter space as discussed by five points in Section 4. 
The darker regions indicate a higher /nl- The red dashed contours are for /nl = 1, 5, 10, 20 from the 
top right to the bottom left. 



are the mass ratio m^/m^ and N c . For this model the energy density ratio at the phase 
boundary is given by 

r = 41 (4.12) 

so we will use r instead of m^/m^ for the parameter space, together with iV c . 

After transforming the parameter space, the free parameters reduce from m^, m a , cro, N c 
to r and N c . The other two degrees of freedom are absorbed in which is fixed by obser- 
vation, and do which does not change the cosmological predictions. 10 The transformation 

10 In fact cro changes the inflation energy scale and the ratio of the curvature perturbations by the inflaton 
and the spectator, but since for ao ^> 4>c, the spectator totally dominates the curvature perturbations and the 
inflation energy scale is too low to be seen from the primordial gravitational waves, its value does not matter 
any more. 



- 12 - 



in the parameter space has the following relations 



3M*P C 

ml = — p \ , (4.13) 

64(2iV* + 1)7t<7q 

The local bispectrum in Eq. (4.8) now simplifies to (after neglecting the O(l) coefficient) 

/NL ~ (2N c + iy r - (415) 

The two dimensional parameter space for r and N c is constrained by the following 
conditions. 

1. The Hubble exit of the pivot scale, the phase boundary and the end of inflation are all 
well separated. So typically we choose 3 < N c < — 3. 

2. The a field should be subdominant, so r <C 1. 

3. The inflaton should provide suppressed curvature perturbations compared to that of 
the spectator, which means cr 3> (/>*■ 

4. The first order slow roll parameter is smaller than unity when a stays on the flat 
potential, i.e. e CT * < 1 and e CTC < 1. 

5. The quantum fluctuations of a should not dominate over its classical slow roll. This 
means Eq. (3.42) is valid. 

Because of the last constraint, the energy density ratio r can not be too small at the phase 
boundary. Therefore the maximum possible value /n? ~ 20 is achieved when 7V C = 3 and 
r ~ 1CT 3 . 

Under the above conditions, we can calculate the spectral index n s , its running dn s /d\n k, 
the local bispectrum /nl, and the trispectrum g^, according to Eq. (3.30), Eq. (3.34), Eq. 
(4.8) and Eq. (4.9). A specific example of /nl is shown in Fig 2, for the parameters iV* = 50, 
o"o = 10Mp and w = 1/3. The regions violating any of the above five conditions are excluded 
and shown by the yellow shaded region in Fig 2. We can read out from Eq. (4.13) the mass 
for <fi here is ~ 2.1 x 10~ 7 M p , so the curvature perturbations from cj) are indeed negligible. 
In addition, for these parameters we have n s = 0.98, dn s /dlnk = —3 x 10 -4 and g^L ^ 1) 
all of which hardly depend on the choice of iV c or r, and fall within the observational bound. 

From Fig 2, we see that the parameter space is limited. In particular, r is constrained 
on both the sides because we need a to be subdominant and its quantum fluctuations not to 
overcome the slow roll motion. Moreover, a hardly contributes to n s , /nl or $nl because all 
of its slow roll parameters are tiny. With the inflaton <j) being the only contribution, we get 
n s ~ 0.98, and small running, bispectrum and trispectrum. For these parameters, we can see 
the local bispectrum strength /nl indeed has a maximum value around 20. The parameter 
space for /nl > 10 is very limited. In this case the major contribution to /nl comes from the 
conversion of the entropy to the curvature perturbations, which becomes non-Gaussian after 
the spectator ends slow roll and decays into a perfect fluid, even though this non-Gaussian 
conversion only lasts for one e-fold or so before the entropy perturbations are redshifted away. 
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5 An inflection point potential for the spectator field 



5.1 Generic inflaton potential 

Flat directions naturally arise in string theory or supersymmetric theories [5]. These flat 
directions can also be candidates for U(a). In most cases, such flat directions can be written 
locally as an effective scalar potential in the form [20, 21] 



where 



Aa = a-a . (5.2) 



Therefore at Act = 0, i.e. a = cro, we will get an inflection/saddle point where U = Uq and 
U" = 0, and Act is the deviation from the inflection point <to- For inflection and saddle points 
we will have 71 > and 71 = respectively, and we always have 73 > 0. In general the higher 
order terms in the effective potential Eq. (5.1), e.g. Act 4 , also provide a small contribution 
to the potential or its derivatives. Here we assume their contribution vanishes for the sake 
of simplicity. 

In this respect the motion of a can be solved as follows. We first obtain the deviation at 
the phase boundary Acr c , from the breakdown of second order slow roll condition r] ac = — 1 

8irM p (U c + V c ) 

Acr c = . (5.3) 

I3U0 

With this we can introduce a very helpful parameter 70 , which tells us how "flat" the potential 
is at the inflection point 



= Hi 2M I [jTft Uo ,,,, 
70 ~ V M p 7sAa c 2 V 2 4n(U c + V C )' 1 ' ' 

Therefore the ratio of U'(a) between the inflection point and the phase boundary is 70/ (1 + 
70). For the inflection point potential we discuss here, typically 70 *C 1. 

As long as we specify the inflaton potential and N c , we are able to solve the equation 
of motion for Act. The slow roll approximation in phase I gives the l.h.s of Eq. (3.3) as 

Aa * dAa M 2 \~^T~ Hyjf Act 

arctan 



Aa c U' U V 7173 V 2 7i M P 

Ml ( 1 



47T7o(t/ c + y c 



(arctan h arctan— ), (5-5) 

V 7o 70/ 



where x = Act/|Act c | is the relative displacement from the inflection point. 

Because a always subdominates the energy density, we can neglect its contribution to 
the Hubble rate and solve the background evolution for 4>{N). By equating Eq. (5.5) (as the 
l.h.s of Eq. (3.3)) with the r.h.s of Eq. (3.3), we derive the evolution of cr by x*(N*) 

x* . 1 ^lo{Uc + Vc) [ m # R v 
arctan — = — arctan 1 — -= / — . (5-t>) 

7o 7o J dr V 
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We can also derive the slow roll parameters that are needed to calculate the cosmological 
observables 



while all the higher order slow roll parameters vanish. When 73 is large and 70 <C 1, we have 
s a ~ 1, and the final results are simplified to 

n s — 1 = —260* + 2x*, (5.8) 



cilnA; 2 



;/„ - l) 2 + 2e^(2r]^ - 3e«^) + 7o - a;*, (5.9) 



10(5 + 6^)6^(26^-^) 5x* 

/nl ~ f1 77— s ^— , (5. 10) 

9(1 + ty)r or 

50(5 + 6w)e^ 25(3x 2 - 7o 2 ) 
5NL 81(l + u;)r2 + 108r2 " (fU1) 

Typically, we have ~ O(10~ 2 ), and e^ c ~ r\$ c < 10 _1 . For this reason in Eq. (5.10), 
the x* term dominates in general, and becomes the major contribution to /nl- The x* term 
comes from r) a * in the original equation Eq. (3.36). A brief estimation for the maximum local 
bispectrum can be achieved for |x*| ~ 10~ 2 with a small r ~ 10~ 4 , so |/nl I ~ 100. 

5.2 Quadratic slow roll infiaton 

As a simple example, we consider the infiaton potential to be 

V(4>) = m 2 (j) 2 , (5.12) 

and we assume it dominates over the spectator potential U(a). After taking the pivot scale 
as iV* = 50, we have a total of 5 free parameters m, Uq, 71,73 and N c . By using the same trick 
as in Section 5, we can fix the overall energy scale to match P^, and switch the parameters 
to r, 70, 73, N c . From Eq. (5.8) to Eq. (5.11) and Eq. (5.6), we can see both the cosmological 
observables and the solution x*(A^) are independent of 73 when 73 is large. Therefore, here 
we take a large 73 and further reduce the parameter space to r, 70, N c . 

The background evolution of <fi can be worked out as a function of ./V 



<XA0 = M p y^±^. (5.13) 
This reduces Eq. (5.6) to 

arctan * = - arctan 1 + 1^±H ln ^±1. (5 .14) 
70 70 4 2JV e + l 

In Eq. (5.14) when iV* = 50 and 70 are fixed, the relative displacement x* is maximized 
when the last term in Eq. (5.14) reaches maximum at N c ~ 18. Since the field a contributes 
to the spectral index by 2x* as in Eq. (5.8), the spectral index would also reach maximum 
at N c ~ 18. The local bispectrum has a — 5x*/6r contribution, so it should reach minimum 
at N c « 18. For the m 2 (j) 2 infiaton potential we have ^ = 0, so <?nl would be positive for 
x 1 > 7o/3 an ^ negative for x 2 < 7q/3. 
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Figure 3. Cosmological observables and 71 are drawn for the inflection point spectator field with 
a quadratic inflaton potential. The x axis is the logarithmic of r, the energy density ratio of the 
spectator field w.r.t the total at the phase boundary. The y axis is the number of e-folds of inflation 
from the phase boundary "c" to the end of inflation. The shaded yellow regions are excluded by 
the five constraints discussed in Section 4. The shaded green regions are observationally favored by 
spectral index within n s = 0.9608 ± 0.0160, whereas the solid green lines indicate the central value. 
[1] 
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To verify our analytical analysis and check the observational bound for the parameter 
space, we take the specific case of TV* = 50, 73 = 10 10 , 70 = 0.15 and w = 1/3. As the 
power spectrum can be determined by the overall energy scale, we plot the rest of the CMB 
observables w.r.t the r and N c coordinates in Fig 3. The parameter 71 is also shown in the 
transformed parameter space. The parameter space has the same exclusion conditions as 
discussed in Section 4 which are shaded yellow. Here the energy scale for the spectator field, 
U Q , varies a lot from 10" 26 M^ to 10" 17 M^. 

In Fig 3, we find agreement with the analytical analysis. The spectral index drops 
out of the "2cr" confidence level around N c ~ 18 because its maximum value exceeds the 
observed central value. Also at N c ~ 18, the local bispectrum reaches its minimum in the 
same time. The running of spectral index is typically small. The local bispectrum and 
trispectrum can both attain large or small values with plenty of parameter space within the 
current observational bound [1, 35]. In this case, the significant /nl an d <?nl come from the 
rj a */r term in Eq. (3.36). 

6 Conclusion 

We have derived analytically the influence of a generic spectator field a on the temperature 
anisotropy of the CMB. We evaluated the relevant cosmological perturbations seeded by the 
spectator field by 5N formalism, and obtained possible non-Gaussian signatures. Signifi- 
cant local bispectrum and trispectrum on the temperature anisotropy of the CMB can arise 
from the possible curvature of the spectator potential and the dilution process of <r's decay 
products. Since the spectator field decays well before inflation comes to an end, the decay 
products of the spectator field are completely diluted away during inflation, and will not 
generate any isocurvature perturbations. In this respect the spectator field might come from 
the visible sector or from the hidden sector, and need not have any couplings to the Standard 
Model fields 11 . 

We have considered two examples for the spectator potential where it can match the 
current observations. In the first case we have considered a step potential for a and for 
inflaton we have assumed a quadratic potential. In this case /nl ^ 20 and <7nl is very small. 
The model predicts the spectral index to be around 0.98 with a negligible running. In the 
second example we have considered a simple potential for a with an inflection point. The 
/nl and #nl can be as large as the upper bound allowed by the current limit [1], while the 
spectral index and its running are both well within the observational bound. 

To summarize, we have found a new mechanism where the spectator field can indeed 
provide dominant curvature perturbations, with or without a deviation from Gaussian fluc- 
tuations. This suggests that fields that decay during inflation can still play an important 
role, and hence should not be neglected in many plausible scenarios. For any beyond-the- 
Standard-Model physics that introduces plethora of hidden sectors, even though they might 
be dormant and decay during inflation, given sufficient conditions laid down in this paper, 
they can still leave their imprints in the sky. 

11 If the spectator oscillates or its decay products become nonrelativistic before the matter-radiation equiva- 
lence, we need sufficient e-folds of inflation after the spectator ends slow roll, i.e. N c , to prevent any observable 
isocurvature perturbations. In the extreme case where the spectator decays directly into nonrelativistic dark 
matter, we typically require N c > iV»/3 ~ 20 for iV» = 60. 
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A Notations 

A list of notations which are present in the final expressions can be found in Table 1. 



Table 1. Notations 



Symbol 


Definition 


Expression 




Value at the Hubble exit 




c 


Value at the phase boundary 




e 


Value at the end of inflation 




/ 


Phase boundary condition 


M 2 U" + 8vr([7 + V) 


W 


Equation of state during phase II 




r 


Energy density ratio 


U C /(U C + v c ) < i 


s 




(U C + V C )/(U* + V*) < 1 






Uc/U* < 1 


Slow roll parameters (similarly for (ft) 


(-a 


First order 


M*U r2 /16n(U + V) 2 


f)o 


Second order 


M 2 U"/8ir(U + V) 


to 


Third order 


M^U'U'" /{M{U + V)) 2 


K 


Fourth order 


M*U' 2 TJ m /(8tt(U + V)) 3 




Fifth order 


M*U' 3 U""'/(87r(U + V)) 4 



B The expression of a 



From the definition of a, Eq. (3.20), we calculate the integral by changing the variable to 
n((f> c , eft), or just simplified as n, through the slow roll relation 



Then we get 



g -3(l+u>)n( 



MlY 
„ Ty dn. 
8ttV 

n(<t> c ,(t>e) M 2 







8vrU 



-3(l+»)n dn _ 



(B.l) 



(B.2) 



The term 1/V in Eq. (B.2) can be regarded as a function of n and Taylor expanded at the 
phase boundary point, giving 



1 

V 



1 



n dV 



+ 



n 



dV 



V 3 \ dn 



d 2 V 



dn 2 



Vc 



+ 



(B.3) 



For calculations of n-point correlation functions, it is safe to truncate the series expan- 
sion up to (n — l)-th order. We assume the phase boundary c and the end of inflation e are 
well separated by several e-folds so n((ft c ,(ft e ) > 0(3). In this case we can safely extend the 
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upper bound of the integral to infinity because of the exponential suppression in n. Then 
the value of a can be worked out order by order 

a = a 0) + a (2) + a (3) + • • • . (B.4) 

For the local bispectrum, we need the first two terms in Eq. (B.4), which are 

a« = —, ^ ^, (B.5) 

3(l+w)(l -r) 2 ' v ; 

and 

a (2) ^ ( ^ (B Q) 

where r]^ = M^V" /8tt(U + V) is the second order slow roll parameter for cf> when both fields 
present. 

For the trispectrum calculation, a needs to be calculated up to third order. Similar 
derivations yield the third order component for a 

and here 

Mp/'Y" 

*+ = 647r2(C/ + y)2 ( B - 8 ) 
is the third order slow roll parameter for cj). 
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